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The motion of a driven planar pendulum with vertically periodically oscillating point of suspension and under the action of an 
additional constant torque is investigated. We study the influence of the torque strength on the transition to chaotic motions 
' of the pendulum using Melnikov's analysis. 
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1 Introduction 



The periodically driven planar pendulum is one of the simplest mechanical systems which can exhibit chaotic behaviour. It 
, is usually modelled as a mathematical pendulum with linear damping and external or parametric excitation |1 2|. Since 
I— I ' the dynamics of various systems of various kind in mechanical 0] El or electrical Gl HI engineering can be described by 
' the same or the same type of differential equations as the one for the driven pendulum it is not a great surprise, that many 
methods invented to calculate and control chaotic vibrations of nonlinear systems have been tested on this system where the 
mathematical parameters appearing in the above mentioned equations can be given a well interpretable physical meaning. 

In this note we focus on the planar pendulum with a vertically periodically moving suspension point with an additional 
constant torque |5 6 1. The corresponding differential equation has the following form: 



ml^cf) + kef) A- m{lg — Ix) smcj) — T, x = asin^t^ (1) 



where m and I denote the point mass and the length of the mathematical pendulum and T is the extra constant torque applied 
at the suspension point, a is the amplitude of the vertical excitation x{t) of the suspension point and the corresponding 
• frequency of excitation. 
' Introducing dimensionless time t = ujt, Eq. [^in dimensionless form is 

o 

+ (1 +7l]'2 cos r2'r) sine/) = io, (2) 



^ where the corresponding constants are given by: uj = \J gjl, il' = il/o;, a — k/{ml'^uj), ^ = a/l, ti^ ~ T^/ {mPoj'^). 
'tI^ Our goal is to investigate the dynamics of Eq. |2]concerning the occuiTence of chaotic vibrations using different approxi- 
mations. 

> 

^ . 2 Melnikov's analysis 

H ■ 

' In the first approximation, we assume that damping, excitation and extra torque are all small terms which we scale by intro- 
ducing a small parameter e to Eq. 121 

(j) + ea(j) + (1 + ejV,''^ cos fl'r) sin (f>^eio, (3) 

where ae = a,je = 7, toe = tg, respectively. This allows us to apply Melnikov's method [TjjSJ for the unperturbed symmetric 
potential [2J. The critical value of the excitation amplitude 7c which is given by a simple zero of Melnikov function f7l|8l: 



7c 



, ton 



sinh(r!'|). (4) 



is (for to — 0) very similar to that obtained in [2J, with 7^, differing by a constant coefficient depending on the torque to'. 
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Fig. 1 (a) Dependence of critical value of the excitation amplitude 7c versus fi' for a = 0.1 and three values of a constant 
torque to = 0.0, 0.05 and 0.1 in the case of heteroclinic transition, (b) the same dependences for to — 0.2 using two different 
approximations with heteroclinic - T and homoclinic - '2' transitions. Below the curves the system is regular while above 
transit to chaotic oscillations occurs. Note, different scales in Figs, la and b. 
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Fig. 2 Phase trajectories (where v — (fi) for to =0.2 and two different approximations, with heteroclinic (a) and homoclinic 
(b) orbits (marked by 'het' and 'hom'), respectively. 

In Fig. [It 7c versus fl' is plotted, showing the influence of the torque on the existence of chaotic dynamics for (Eq. |3}. 
Note that these results are obtained for small torque compared to damping to/nS ^ a. 

In the second approximation we treat the case of larger torque that is to/(7r8) ~ a. Here we use the recently developed 
method for non-symmetric potentials |9|. This approach deals with the equation 

+ ea^ + (1 + £7^^'^ cos rj'r) sin - to = + ea0 + [l - to + - 0^/6] + e^yfl''^ cos n'r sin = 0, (6) 

where contrary to the first approximation to is not taken to perturbations terms (contrary to Eq. |3} but included into the 
effective unperturbed potential. Consequently it breakes the potential symmetry and instead of a heteroclinc transition with 
two saddle points, a homoclinic transition with a single saddle point is obtained [5] is obtained (see Fig. 2). 

In Fig.^ we compare results of both approximations. Curve '1' corresponds to the small torque approach with heterocUnc 
transition while '2' represents the results with to included in the effective potential (a homoclinic transition). Note, for 
small frequencies both curves match quite well but for higher frequencies there is a growing discrepancy between the two 
approximations. 
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